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A Polarimetric Extension of the van Cittert-Zemike 
Theorem for Use with Microwave Interferometers 


J.R. Piepmeier, Member , ; 


Abstract — The van Cittert-Zernike theorem describes the 
Fourier- transform relationship between an extended source and 
its visibility function. Developments in classical optics texts use 
scalar field formulations for the theorem. Here, we develop 
a polarimetric extension to the van Cittert-Zernike theorem 
with applications to passive microwave Earth remote sensing. 
The development provides insight into the mechanics of two- 
dimensional inteferometric imaging, particularly the effects of 
polarization basis differences between the scene and the observer. 

Index Terms — Interferometer, polarization, polarimetric, 
STAR 


I. INTRODUCTION 

T HE desire to obtain fine spatial and temporal resolution 
passive microwave imagery of Earth from space drives 
the need for microwave instruments with large apertures. 
As aperture size increases beyond current practice, however, 
imaging by mechanical scanning becomes challenging because 
of structural dynamics. Electronic scanning methods, on the 
other hand, alleviate the need to move large antenna structures. 
One such electronic method is microwave interferometry, also 
called synthetic thinned array radiometry (STAR) [1], An 
interferometer measures the complex spatial correlation field, 
or visibility function, of the incident radiation originating 
from a source. The van Cittert-Zemike theorem describes the 
Fourier transform relationship between a source distribution 
and its associated visibility function. In classical optics texts, 
the van Cittert-Zemike theorem is developed using scalar fields 
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demonstrated STAR instruments have been single polarization 
and can be described by the scalar theory (see e.g., [4]). 
An interesting effect occurs in two-dimensional (2-D) STAR 
imaging: the observer and source polarization basis coordi- 
nates do not match in all directions of observation. Thus, the 
polarization basis can be rotated differently for each pixel in 
an image. This indeed is recognized by the developers of the 
2-D STAR instrument on the Soil Moisture Ocean Salinity 
(SMOS) mission [5]. Here, the polarimetric extension to the 
van Cittert-Zemike theorem lends particular insight into this 
process. 

Our development follows the form of [2] with the addition 
of partially polarized radiation from the extended source. 
Another recent development included partially polarized light 
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Fig. 1. Infmtesimal source m within the extended source a radiates and is 
observed at the far-zone location R . The observer polarization basis is rotated 
by a with respect to the source. 


[6]; however, the analysis used the paraxial approximation 
for propagation in beam optics. Here, as in [2], we apply 
the spherical wave approximation appropriate for far-field 
propagation of waves originating from infinitesimal point 
sources. The polarimetric visibility function is expressed using 
the so-called coherency vector [7], which contains the same 
information as the mutual coherency matrix used in [5] and 
the beam coherence-polarization matrix used in [6]. The 
coherency vector can be transformed to the modified Stokes 
visibility vector, common in radio-astronomical polarimetry 
[8], by a simple matrix operator. Additionally, we consider 
the differing polarization basis coordinates of the source and 
the observer. Polarization rotation of the electric fields is 
performed using Jones matricies. (For a formal description 
relating Jones matricies to the Stokes vector see [9].) Our result 
is instructive for interpreting 2-D STAR imagery of the Earth. 

II. Radiated Field 

The far-zone electric field E m at location R arising from 
an infinitesimal element m within the extended source a can 
be written (see Fig. 1): 
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where R = | R | is the range, k is the direction of propa- 
0%^jaop4clSEEhe speed of light, and / is the frequency. This 
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expression is the vector phasor representation of a partially- 
polarized quasi-monochromatic field with complex envelope 
A m . The complex vector random process A m contains both 
the direction and polarization dependence of the radiation 
emitted by the surface element and is assumed to be zero mean 
and stationary. The electric field is related to the brightness 
distribution due to m by [10] 
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where 7] 0 is the intrinsic impedance of free space and dfl is a 
differential solid angle. Its envelope is similarly related by 

<3) 


Br, 


E m = 


rWn,p 

Em,q 


Am 


A 


rn.jp 

J 


(k) da = 


Ip 

h 

U 


da = 


Vo 


J^m,p 

01 


01 


2 Re < 
21m < 


Amjp f k j 




m,q I 


m,q 1 


P(a) = 

such that 

X 


’ P P' 

q p' ' 


COS Ot 

sin a 

p-q' 

q? . 


— sin a 

cos a 


= P (a) E„ 


A' 


— P ( a ) A v 


where a is the angle of polarization basis rotation (PBR): 
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III. Visibility Function 

The coherency vector of electric fields received by two 
isotropes with displacement D is defined by [7]: 

7(d) = (X (■£,*) ®X* (r + D,t)) (10) 

where ® is the outer product operator. Use of the coherency 
vector to represent the visibility function is convenient because 
its elements relate directly to the co- polarized and cross- 
polarized complex correlations measured with an interferom- 
eter. The outer product is expanded to yield 


where da — dPt/R 2 is differential surface area. Note, the 
above expression is the radiation intensity or pattern of the 
infinitesimal source. 

The electric field and its envelope have two polarization 
components in directions p and q, the polarization basis 
vectors of the source frame: 
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The total electric field received by an isotrope is found by 
summing over all elements within the extended source: 
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Substituting (9) into (10) and using the parallel-ray approxi- 
mation, we obtain 
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The polarization basis vectors form a right-handed Cartesian 
basis with k such that k = p x q. The process A m is related 
to the modified Stokes vector I a ^kj of the source in the p-q 
polarization basis: 
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A property from tensor algebra allows some rearranging: 
7 (d) = £ (13) 
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where the subscript a is used to denote the continuous 
dependence upon location within the source [10]. 

The polarization basis vectors in the observation frame p' 
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to the source frame by a Jones matrix of polarization basis 
rotation: 


(Am (k m ) ® A* m (k m ) ) 


R R + D 


The time differences in the expected value are neglected 
because A m is assumed to remain self-coherent at a time lag 
of D/c. Examining the expected value, we hnd 


(■ A m (k m ) ® Kn (km) ) = VoSI a (k m ) da (14) 


(6) where 
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transforms the modified Stokes vector (5) to a coherency 
vector representation [7]. The coherency vector can now be 
rewritten as an integral over the source 

(16) 

or over solid angle 
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The outer product P (a) ® P (a) acts as a PBR matrix in 
coherency vector space. Equation 17 can be written with the 
familiar PBR matrix for the modified Stokes vector [1 1] 


V ( d) SR (a) I, (k) e^-^dn 

where 

cos 2 a sin 2 a -sin 2a 0 

sin 2 a cos 2 a — | sin 2a 0 
— sin 2a sin 2a cos 2a 0 
0 0 0 1 



(18) 


(19) 


The expression in (18) is the polarimetric extension of the 
classical van Cittert-Zemike theorem. 


IV. Application to Earth Remote Sensing 


We can consider the visibility fimction and the brightness 
distribution to be a Fourier-transform pair described by (18): 

V (d) - rjoTt (k) (20) 


where 


t (k) = R (a) I a (k) 


( 21 ) 


is the brightness distribution in the observer polarization 
basis. Thus, the key difference between the polarimetric and 
scalar formulations of the van Cittert-Zemike theorem, besides 
(18) being in vector form, is the polarization basis rotation 
between the source and observer. This variation impacts the 
interpretation of 2-D STAR imageiy of the Earth. For example, 
when the measured visibility function is inverted to obtain 
an image, each pixel in the reconstructed image will have 
a polarization basis different from the Earth. To obtain the 
brightness distribution in the Earth’s polarization basis, the 
inverse of the polarization basis rotation matrix is applied for 
each observation direction: 


I a (k)=R \a)T a (k). (22) 

Here, we consider the case of a nadir-viewing 2-D STAR 
and will determine the polarization basis rotation angle across 
the image. We define the Earth’snatural polarization basis such 
umi die hurizoniai poiarization h e is parallel to the surface and 
the vertical polarization is v e = h c x k with k pointing away 
from the surface. As shown in Fig.2, the origin of the observer 
coordinate frame is placed at height H above the surface with 
z pointing down and k = — r.The source polarization basis is 
aligned with Earth’s natural polarization basis, which yields 

p = v e = 0 q = h e = - (23) 


The observer polarization basis is arbitrarily aligned to obtain 
p' = x at nadir (6 = <j> = 0). The conventional Ludwig’s third 
definition for polarization basis of an antenna is used to define 
[ 12 ] 

p' = cos <j> 0 — sin <j> <j> (24) 

cf = — sin cf) 0 — cos <j> (f) 

Finally, the polarization basis rotation angle for a pixel in 
direction (6, <p) is found by substituting (23-24) into (8): 

a = 0. (25) 



Fig. 2. With the observer coordinate frame located at H above the Earth, the 
polarization basis is rotated by a for pixel m located in the (i 9 , <j>) direction. 


Thus, the poiarization basis in a single image channel will 
rotate with the azimuth angle in the image. For example, the 
/p image will contain a mixture of I v , I h , and U depending 
upon the azimuth angle. With <j> = 0, I' p = I V9 but with 
</> = 90°, ij, = In- Thus, to obtain the Ih image, for example, 
the entire I ' vector is needed. 

v. Discussion 

Here, we developed the polarimetric version of the van 
Cittert-Zemike theorem, which describes the Fourier transform 
relationship of a brightness distribution with its visibility 
function. Using the coherency vector to represent the polari- 
metric visibility function yields a compact expression that pro- 
vides insight into interferometer operation. For example, this 
expression demonstrates how polarization basis differences 
between source and observer affect a reconstructed image. 
When imagery is desired in the Earth’s natural polarization 
basis, the full complex mherency vector needs to be measured. 
Designers, however, may be tempted not to measure cross- 
polarization correlations to save engineering resources. If only 
the two co-polarized visibilities are measured, then PBR can- 
not be completely removed. Mathematically, R becomes a 2x2 
matrix that is ill-conditioned near ±45°, at which it becomes 
singular. On the other hand, imagery in the instrument’s 
polarization basis may be acceptable. In atmospheric sound- 
ing, legacy real-aperture instruments already exhibit varying 
polarization basis. Thus, instrumental PBR, as described by 
the polarimetric van Cittert-Zemike theorem, could be in- 
corporated into the retrieval of atmospheric parameters from 
interferometer imagery. 
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